Driven Tunneling Dynamics: Bloch-Redfield Theory versus Path Integral Approach 
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In the regime of weak bath couphng and low temperature we demonstrate numerically for the 
spin-boson dynamics the equivalence between two widely used but seemingly different roads of 
approximation, namely the path integral approach and the Bloch-Redfield theory. The excellent 
agreement between these two methods is corroborated by a novel efficient analytical high-frequency 
approach: it well approximates the decay of quantum coherence via a series of damped coherent 
oscillations. Moreover, a suitably tuned control field can selectively enhance or suppress quantum 
coherence. 
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The dynamics of driven quantum systems which inter- 
act with a large number of environmental degrees of free- 
dom ||lj-|3|] plays an increasingly prominent role: its vast 
appHcability ranges from tunneUng phenomena in solid 
state physics, the study of electron and proton transfer 
in condensed phases, to the gate operation in quantum 
computing devices H, to name but a few. In particu- 
lar, the use of properly tailored external driving forces 
enables one to selectively manipulate a quantum trans- 
port process. The various communities typically rely on 
different methods of description. The two most popu- 
lar approaches for a portrayal of the time evolution of 
the corresponding reduced density matrix (RDM) are ei- 
ther based on the system-bath coupling expansion ob- 
tained by use of a projector operator method - being 
commonly known under the label of the Bloch-Redfield 
formalism ~ or on the expansion in the coupling matrix 
element A (such as a tunnel splitting) by use of (real- 
time) path integral methods. Nevertheless, there exists 
practically little crosstalk between the practitioners of 
the two approaches, and even more, not much of detailed 
comparison between the two seemingly different roads of 
approximation needed for practical calculations. 

For the archetype quantum system of a driven spin- 
boson dynamics, namely the driven dissipative two-state 
system (TSS) dynamics 0, the application of the so 
termed noninteracting blip approximation (NIB A), i.e. 
the leading order result in the tunnel coupling A^, pro- 
duced many impressive successes in entangling the com- 
plexity of driven open quantum systems. This scheme 
works best in the regime of strong friction and/or high 
thermal temperatures. Much less is presently known, 
however, about the corresponding complexity of the 
driven dynamics in the deep quantum regime at low 
temperatures and weak system-bath coupling, where the 
NIBA is failing and higher order terms in the series in 
A must be accounted for [g|,||. In practice, this latter 



regime is of relevance for many situations such as, e.g., 
for the challenge of "battling decoherence" in quantum 
computing schemes Q. 

Our main objective with this work is to enlight the 
advantages and disadvantages of the two approaches. In 
doing so we present three major findings: (i) We numer- 
ically demonstrate the equivalence for the driven tunnel- 
ing dynamics between the path integral method beyond 
NIBA and the coupled set of nonstationary, Markovian 
Bloch-Redfield equations, (ii) Starting from the gener- 
alized master equation (GME) for the RDM, obtained 
within the path integral approach, we arrive at a novel 
analytic high driving frequency approximation that com- 
pares well with comprehensive numerical findings, (iii) 
With this analytical result one can efficiently determine 
the optimal control of quantum coherence. 

Our starting point is the driven spin-boson Hamilto- 
nian 0] where the TSS is bilinearly coupled to an ensem- 
ble of harmonic oscillators, i.e.. 



H{t) - -h[Aa, + €{t)a,]/2 
+ ^;iw,(6j6, + l/2)4 



<t.^q(6, + S1)/2, (1) 



with cTi being Pauli spin matrices. Here A describes the 
coupling between the two states, and e(i) is the exter- 
nal, time-dependent control field. The basis states are 
chosen such that \R) (right) and \L) (left) are the lo- 
calized eigenstates of the "position" operator ct^. All 
effects of the Gaussian bath on the TSS are captured 
by the force autocorrelation function ||l|-^ M{t) = 
- j„ divJiuj) ,i;h(L/2fc^T) ^ ^here the spectral den- 
sity of the environment, J{uj) — ttH^ ^,- cfS{oj — uji) = 
2TTauje^'^^'^'', is assumed to be of Ohmic form with ex- 
ponential cutoff and dimensionless coupling strength a. 
The dynamical quantities of interest are the expectation 
values fTi(i) '.— Tr{/5(i)iTi} which, together with the unit 
matrix /, comprise the complete reduced density matrix 



p[i) = 1/2 + Y,i=x,y.z'^ii^)^il'^- In the following we as- 
sume that at time i = the particle is held at the right 
site fjz — +1, with the bath being in thermal equilibrium. 
Path-Integral Approach. ~ For a harmonic bath the 
exact formal solution for the evolution of the iJi{t) can 
be expressed in terms of real-time double path integrals 
[0-p|. This procedure yields the formally exact set of 
equations HJI,^ 



^z{t) 



(^x{t) 



dt' 



dt' 



Ki-\t,t') ~ Ki+\t,t') a^it') 
Ki+\t,t')+Ki-\t,t')azit' 



(2) 



and <7y(t) = —&z{t)/lS.. Here, the kernels Kl ', i = x,z 
are found in the form of a series expansion in A. Be- 
cause the exact series expression cannot be evaluated to 
all orders, approximation schemes necessarily must be in- 
voked. A familiar scheme is the noninteracting-blip ap- 
proximation (NIB A) f^-y], which corresponds to a trun- 
cation of the series expansion to lowest order in A. The 
NIBA is approximatively valid only for the dynamics of 
(Tz{t) if on average (e(<)) = 0. However, in the presence 
of a static asymmetry component, it breaks down for 
weak damping and low temperatures P,p[ . A systematic 
weak damping approximation for the kernels K^ in (y), 
which circumvents the weaknesses of the NIBA has been 
discussed in |^,|[ . By keeping track of the bath- induced 
correlations to linear order in a, the whole series expan- 
sion in A can be summed analytically. The kernels Kz 
read 

Ki+'^ {t, i') = A^ cos[C(i, t')] [1 ~Q\t~- 1')] 

+ [ dt2[ dtiAUm[C{t,t2)]Poit2,h)sm[Citi,t')] 
X [Q'{t - t') + Q'{t2 - ti) - Q\t2 - t') - Q'{t - ti)] , 

ifi-)(t,i')-A2sin[C(i,t')]Q"(t-0 

t /•*2 

dt2 dti A4sin[C(t,t2)]/'o(i2,ti)cos[C(ii,t')] 
t' Jf 



x[Q"{t-t')-Q"{t2-t')]. 



(3) 



Here, the first term in Kz represents the weak-coupling 
form of NIBA. In the remaining contribution the term 
Pa{t2,ti) accounts for all tunneling events during the 
time interval [^1,^2] that are not influenced by damp- 
ing. Hence -Po(i2,^i) solves the generalized master equa- 
tion (GME) for crz(i) (0) with the zero-damping kernels 
K^+){t^t') = A2cos[C(i,t')] and K^'\t,t') ^ 0, where 
C{t,t') — Jj, dt"t{t") captures the effects of the external 
force. The bath-influence is encapsulated in the functions 
Q'{t) and Q"{t) being the real and imaginary part, re- 
spectively, of the twice integrated bath correlation func- 
tion 7W(i). 



Bloch - Redfield Formalism. - In the Nakajima- 
Zwanzig theory |7| it is well known how to construct an 
exact generalized master equation for the reduced den- 
sity matrix with the help of projection operators. For 
intermediate to high temperatures and/or strong damp- 
ing, but for arbitrary driving, a master equation for p{t) 
can be obtained within the small polaron theory, yield- 
ing equations which are equivalent to the NIBA MM . For 
weak coupling to the bath the projection operator tech- 
nique yields the GME in Born approximation that can 
be further simplified to the Markovian kinetic equations 
without loss of accuracy to the leading order in dissipa- 
tive coupling. For strong harmonic driving this objective 
was first achieved in 1964 by Argyres and Kelley [O. 
Following the reasoning in ||ll| the kinetic equations for 
the RDM of a stochastically driven TSS were found in 
[12a] and in a different way in [12b]. Generalizing [|ll|jl2| 
to the case of a spin-boson problem with an arbitrary 
control field we find the coupled equations 

&xit) = e{t)(7y - T^x(t)ax - Ta;zit)az - Ax{t) , (4) 

^yi^) = -e(i)o-x + Actz - Tyy{t)ay - Tyz{t)az - Ay{t) , 

with Tyy{t) — Txx{t) and az = — Actj,. Here the 
time-dependent rates Tij(t) = J^ dt'M'{t — t')bij{t,t'), 
together with the inhomogeneities Ax{t) ~ ImF(i), 
Ay{t) = ReF{t), with F{t) = 2j*dt'M"{t - 
t')URR{t,t')URLit,t') determine the dissipative action of 
the thermal bath on the TSS. The functions M' and 
Ai" are the real part and imaginary part, respectively, of 
the correlation function A4. The quantities URR{t,t') — 
{R\U(t,t')\R) and URL{t,t') = {R\U(t,t')\L) are ma- 
trix elements of the time evolution operator U{t,t') of 
the non- dissipative driven TSS. The functions 6^ read 
bxx = \Urb,\^ - IUrlI"^, b^z = 2 RgUrrUrl, and byz = 
—2 IuyUrrUrl. This main result in Eq. (^ yields for the 
first time a consistent Bloch-Redfield-type description of 
the externally driven spin-boson dynamics. Equations of 
the form (H) were derived by Bloch and Redfield in 1957 
|13| to describe spin relaxation in nuclear magnetic reso- 
nance, and in [|14| for the dynamics of the undriven spin- 
boson problem. Our set of Eqs. (Q) generalizes |1J] to 
general driving forces. Note that these derived equations 
are valid in the parameter region aln(a;c/A) <C 1, where 
the frequency corrections to the dynamics incurred due 
to the dissipation are small, and the perturbative treat- 
ment is fair. One can show that for the undriven case, 
e{t) — eo, the analytic solution of Eq. (Q) in first order in 
a reproduces the analytical path integral weak-damping 
results, cf. [g,|| and (||) below with zero ac-field. 

Analytic high frequency solution. - Up to here no as- 
sumptions on the deterministic control field have been 
made. Next, we focus our attention on a monochromatic 
field of the form e{t) ^ ea + s cos Oi. Moreover, we restrict 
our investigations on the CTz(i)-dynamics, as this quantity 
is of prime interest for describing tunneling properties. 




FIG. 1. Matching between Path Integral and Bloch-Redfield. 
The comparison of the dynamical Eqs. (H), (0) and (Bl) for 
unbiased TSS-dynamics depicts excellent agreement. For this 
resonant (eo = n*Q.,n* = 0) condition the dynamics is well 
described by (H) with the single- mode frequency Oq — Aq. 
Here and in the following figures frequencies are expressed in 
units of A, times in units of A^^. The temperature is zero 
throughout. 

Because the path integral approach yields a closed 
integro-differential equation for az{t), we start from the 
generalized master equation (0). In the high frequency 
regime {Vt :$> {A,eo,rfl}, with r^ defined in (^ below) 
a good approximation to the dynamics of az (t) amounts 
to perform the substitution 7^[C(i,i')] — ' (^[C(i,*')]) = 
Jo{x)T^[f^o(t - 1')] into the kernels Ki'^\t, t'), where 7^ = 
cos or sin and x = |f sin ^*~* ' . Here ( ) denotes time- 
averaging, and Jo is the zero order Bessel function. The 
resulting generalized master equation is in the form of 
a time-convolution. A solution is conveniently obtained 
by use of Laplace transformation techniques, upon gen- 
eralizing a line of reasoning proposed in Ref. ||5|]. By 
expanding the Bessel function Jo(x) in Fourier series, 
and upon introducing the field-dressed tunneling split- 
tings A„ — |J„(s/57)|A and the photon- induced asym- 
metries e„ = eg — 7T,J7, we end up with 



^z{t) = Poc + {Vo - Po.)e-^''' +YCn cos{9„t)e 



(5) 



Conservation of probability yields Vo + X^n^" ~ 1' ^^^^ 



m^n 



The damping rates and the averaged nonequilibrium 
value Pnn read llisi 



r_R — 2 2^r„ , r„ — -c„ /„ s{6n) 

n 



A2 



(6) 

(7) 



Here /„ = V^ ^« E™ Am/[£™(^' " 4)], and S{9) = 
J{9) coih{'h6 / 2k bT) . The infinite set of frequencies 0„ is 
determined by the pole equation for the undamped TSS 



n(^^"^')+E^'n(^™-^')=o 



(8) 
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FIG. 2. Driving induced quantum coherence phenomena. In 
the presence of a quasi-resonant high frequency field away 
from the zeros of J„*(s/r2), the population difi'erence cyz{t) 
exhibits a coherent oscillatory decay which is dominated by 
a single mode oscillation frequency On* . A comparison be- 
tween the predictions of the analytical solution (0) with just 
the single mode frequency 9\, for a near-resonant field (i.e., 
n* = 1 with ei = |eo - n| = 0.2A) with the Bloch-Redfield 
result in M) is depicted. Note that in the undriven situation 
(s = 0) the TSS dynamics is almost completely localized. 

Finally, to approximately take into account bath- 
induced frequency-shifts the tunneling frequencies 0„ are 
evaluated from (0) upon substituting A„ — > A„(l — 
a\n(ujc/ ^)) '■= A„. Thus, in this high frequency regime 
the system generally still exhibits damped coherent oscil- 
lations, as in the undriven case, although, an infinite set 
of oscillation frequencies 0„ with corresponding damp- 
ing rates r„ enters this driven dynamics. Superimposed 
to these coherent oscillations there occurs an incoherent 
decay with rate F^ towards Poo • 

In figures 1-3 we depict comparisons amongst the nu- 
merical predictions of the Born-Markov equations (|j), 
the path-integral GME (g), and the analytical solution 
(H) for small Ohmic friction and zero temperature. For 
the driven dynamics the agreement is remarkable. It in- 
creases further with increasing temperature (not shown) . 
To achieve a convergence of (pT) a truncation of the pole 
equation (g) to five (or less, cf. Figs. 1,2) modes, charac- 
terized by (A„,e„), turned out to be sufficient. In Fig. 
1 the infiuence of an unbiased (eg = 0) control field is 
investigated. This corresponds to a resonant (eq = n*Q) 
field with n* = 0. In Figs. 2 and 3 the case of a fi- 
nite bias eo 7^ is depicted. Fig. 2 depicts the near- 
resonant situation Jeq — n*D,\ — |e„. | ^ {A„», |eo|} away 
from the zeros of Jn'{s/^): the coherent dynamics is 
now already well captured by the single resonant mode 

frequency 0„* ~ y A^j. -I- e^^. . This finding generahzes 
the small-dc-bias analysis in |10[| . In addition, we de- 
duce from the parameters chosen in Fig. 2 that our ap- 
proach can even work for intermediate driving frequencies 
(il w eq). Due to the fact that A„. < A, and e„* < eg, 
the field-induced oscillation frequency 0„* can be much 
smaller than in the undriven case. In the off-resonance 
situation (|e„| > A„ for all n) of Figs. 3a and 3b, <7z{t) 
exhibits a complex interference pattern with 
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FIG. 3. Controlling tunneling. In the presence of an 
off-resonance no net separation of time scales occurs and the 
population (Tz(i) shows a complex interference pattern (Fig. 
3a). - Note that the numerical solutions of Bloch-Redfield 
and path-integral equations coincide within linewidth. The 
TSS dynamics is dominated by an incoherent decay towards 
its asymptotic limit (Fig. 3b), so that quantum coherence is 
lost. The incoherent decay rate F/j, however, can be strongly 
diminished. This is demonstrated in the upper left inset where 
the photon assisted decay rate Tr is plotted vs. the dc-bias eo. 
It exhibits characteristic resonance peaks at multiple integers 
of the driving frequency Q. These peaks are shifted replicas 
of the dc-driven (s — 0) rate with different weights. Thus, 
a suitable chosen bias can enhance or suppress the decay of 
populations. Finally, the lower right inset shows the averaged 
nonequilibrium population difference Poo. It exhibits a non- 
monotonic dependence on the dc-bias when combined with a 
high frequency field. For appropriate values of the dc-field a 
population inversion {Poo < when eo > 0, and vice versa) 
can occur. 

quantum coherence suppressed. Moreover, the decay 
towards the nonstationary equihbriuni value occurs on 
a much longer time scale as compared to the case with 
s = 0. This result, observed recently in iQ, can be un- 
derstood via close inspection of the upper left inset in 
Fig. 3b, where the averaged decay rate Tn is plotted vs. 
eo- For the chosen parameters the decay rate is strongly 
diminished. Moreover, the lower right inset depicts the 
averaged nonequilibrium value Poo vs. the dc-bias eq. 
Here, photon assisted tunneling rules the possible inver- 
sion of asymptotic population (i.e. Poo < 0, for eo > 0, 
and vice versa). 

In concluding we maintain that in the perturbative 



regime {a\ii{LUc/A) <^ 1) it is numerically advantageous 
to evaluate the weak coupling tunneling dynamics by 
use of the nonstationary Markovian Bloch-Rcdficld equa- 
tions, as compared to the non-Mar kovian path integral 
GME. We find numerically perfect agreement as depicted 
with Fig. 1-3. Within the time scale of tunneling we 
find no observable non-Markovian effects. On physical 
grounds, the same remarks apply to the time evolution 
of the full density matrix. Note, however, that for the 
scaling regime (i.e. cutoff cUc -^ oo) Bloch-Redfield the- 
ory increasingly fails; it then is necessary to correct even 
the path integral GME with an additional, concocted 
renormalization scheme B. Finally, our novel analyt- 
ical scheme may prove prominent in order to optimize 

quantum coherence. 
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